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Abstract 

Recent measurements of the properties of cosmic rays above 10 17 eV 
are summarized and implications on our contemporary understanding of 
their origin are discussed. Cosmic rays with energies exceeding 10 20 eV 
have been measured, they are the highest- energy particles in the Uni- 
verse. Particles at highest energies are expected to be only marginally 
deflected by magnetic fields and they should point towards their sources 
on the sky. Recent results of the Pierre Auger Observatory have opened 
a new window to the Universe — astronomy with ultra high-energy par- 
ticles. 



1 Introduction 

The Earth is permanently exposed to a vast flux of high-energy particles 
from outer space. Most of these particles are fully ionized atomic nuclei with 
relativistic energies. The extraterrestrial origin of these particles has been 
demonstrated by V. Hess in 1912 |Hess(1912)| and he named the particles 
"Hohenstrahlung" (high-altitude radiation) or "Ultrastrahlung" (ultra radia- 
tion). In 1925 R. Millikan coined the term "Cosmic Rays". They have a three- 
fold origin. Particles with energies below 100 MeV 1 originate from the Sun 
|Ryan(2005)[|Kahler et al.(2005)| . Cosmic rays in narrower sense are particles 



with energies from the 100 MeV domain up to energies beyond 10 20 eV. Up 
to several 10 GeV the flux of the particles observed is modulated on different 
time scales by the heliospheric magnetic fields [Fichtner(2005) Heber(2005)| . 



Particles with energies below 10 17 to 10 18 eV are usually considered to be of 
galactic origin [Gaisser and Stanev(2006) Gaisser(2006) Strong et al.(2007)| 



Horandel(2008)| |Horandel(2007a)| . The Larmor radius of a particle with en- 
ergy Eis (in units of 10 ib eV) and charge Z in a magnetic field B^q (in fiG) 
is 

r* = 1.08^PC (1) 

yielding a value of tl = 360 pc for a proton with an energy of 10 18 eV in the 
galactic magnetic field {B^q sa 3). This radius is comparable to the thickness 



1 In this review we use the particle physics energy units MeV= 10 6 eV, GeV= 10 9 eV, 
TeV= 10 12 eV, PeV= 10 15 eV, and EeV= 10 18 eV; 1 eV= 1.6 ■ 10" 19 J. 
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of the galactic disc and illustrates that particles (at least with small charge Z) 
at the highest energies can not be magnetically bound to the Galaxy. Hence, 
they are considered to be of extragalactic origin | Na gano and Watson(2006jj 
Bergman and Belz(20"07)||Kampert(2008)| . 



In the present article, we focus on recent results concerning the origin of 
the extragalactic particles. Cosmic rays with energies exceeding 10 20 eV are 
the highest-energy particles in the Universe. Particles at highest energies are 
only marginally deflected in the galactic magnetic fields, following |T]) they 
have a Larmor radius rt > 36 kpc, exceeding the diameter of the Milky 
Way. Thus, they should point back to their sources, enabling astronomical 
observations with charged particles. 

Several questions arise, concerning the origin of highest-energy cosmic 
rays. Among them are: 

- What are the energies of the particles? (Sect.H]) 

- What are these particles? Are they protons, nuclei of heavy atoms like 
oxygen or iron, furthermore are they photons or neutrinos? (Sect. [5]) 

- Where do they come from? Can we learn something by studying their ar- 
rival directions? (Sect.E]) 

- How do they propagate to us? Do they suffer any interactions? (Sect. [4]) 

In the following sections (@] to E]) recent experimental results are com- 
piled and their implications to answer the questions raised above are dis- 
cussed. Before, possible scenarios for the sources of the particles are summa- 
rized (Sect. [2~Ij) and mechanisms are discussed which are important during the 
propagation of the particles through the Universe (Sect. l2.2| . The detections 
methods applied are sketched in Sect. 03 



2 Sources and Propagation 
2.1 Sources 

The energy density contained in the flux of extragalactic cosmic rays can be in- 
ferred from the measured differential energy spectrum dN/dE |Halzen(2006) 

4tt f EdN , 
? E = -Jl3dE dE > (2) 

where (3c is the velocity of particles with energy E. To estimate the en- 
ergy content of the extragalactic component, assumptions have to be made 
about the contribution of galactic cosmic rays at energies in the transition 
region (10 17 — 10 18 eV). The extragalactic component needed according to the 
poly-gonato model |Horandel(2003a)| to sustain the observed all-particle flux 
at highest energies has an energy density of pe — 3.7 • 10~ 7 eV/cm 3 . The 
power required for a population of sources to generate this energy density 
over the Hubble time of 10 10 years is 5.5 • 10 37 erg/(s Mpc 3 ). This leads to 
« 2 • 10 44 erg/s per active galaxy or w 2 • 10 52 erg per cosmological gamma 
ray burst |Gaisser(1997)| . The coincidence between these numbers and the 
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Figure 1: Size and mag- 
netic field strength of possi- 
ble sites of particle accelera- 
tion (Hillas diagram) . Accel- 
eration of cosmic rays up to 
a given energy requires con- 
ditions above the respective 
line |Ostrowski(2002)| . 



observed output in electromagnetic energy of these sources explains why they 
are considered as promising candidates to accelerate highest-energy cosmic 
rays. 

The characteristic size of an accelerating region can be estimated for mod- 
els of gradual acceleration, where the particles make many irregular loops in a 
magnetic field while gaining energy |Hillas(1984)l . The size L of the essential 
part of the accelerating region containing the magnetic field must be grater 
than 2tl. A closer look reveals that a characteristic velocity /3c of scattering 
centers is of virtual importance [Hillas (1984)], which yields the expression 

B^ G L pc > 2E 15 /(Z{3). (3) 

It relates the characteristic size L pc (in pc) and magnetic fields B^q of ob- 
jects being able to accelerate particles to energies Ei$. Several possible ac- 
celeration sites are explored in Fig.[H where the magnetic field strength is 
plotted as function of their typical sizes |Ostrowski(2002)| . The lines ac- 
cording to ([31) represent the conditions for protons and iron nuclei of dif- 
ferent energies, as indicated. Objects capable to accelerate particles above 
a respective energy should be above the respective line. As can be in- 
ferred from the figure, most promising candidates to accelerate highest- 
energy cosmic rays are gamma ray bursts and active galactic nuclei (AGN) 
|Ginzburg and Syrovatskii(1964) Hillas(1984)| . These objects are typically 



in a distance of several tens of Mpc to the Earth. Interactions in the source 
itself or in the vicinity of the source of hadronic particles (protons, nuclei) 
yield neutral and charged pions, which subsequently decay into high-energy 
photons and neutrinos. 
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Alternatively, so called "top-down models" are discussed in the literature 
|Hill and Schramm(1983)l|Olinto(2000)||Bhattacharjee and Sigl(2000j] . They 



have been motivated by events seen by the AGASA experiment above the 
threshold for the GZK effect [Takeda et al.(1998)| . It is proposed that ultra 
high-energy particles (instead of being accelerated, "bottom-up scenario") 
are the decay products of exotic, massive particles originating from high- 
energy processes in the early Universe. Such super-massive particles (with 
mx ^> 10 11 GeV) decay e.g. via W and Z bosons into high-energy protons, 
photons, and neutrinos. 

2.2 Propagation 

On the way from their sources to Earth the particles propagate mostly outside 
galaxies in intergalactic space with very low particle densities. In this envi- 
ronment the most important interactions of cosmic rays occur with photons of 
the 2.7-K microwave background radiation, namely pair production and pion 
photoproduction Hill and Schramm(1985)| . 



On the last part of their way to Earth they propagate through the Galaxy. 
However, since particles at the highest energies travel almost along straight 
lines they accumulate a negligible amount of material during their short travel 
through regions with relatively high densities. Thus, interactions with the 
interstellar material can be neglected. 

The Universe is filled with about 412 photons/cm 3 of the 2.7° K microwave 
background radiation. Shortly after the discovery of the microwave back- 
ground it was proposed that ultra high-energy cosmic rays should interact 
with the photons, leading to a suppression of the observed flux at highest 



energies Greisen(1966) Zatsepin and Kuz'min(1966) |. This effect is called 



after its proposers the Greisen-Zatsepin-Kuz'min (GZK) effect. A nucleon of 
energies exceeding Eqzk ~ 6 • 10 19 eV colliding head-on with a 2.7° K photon 
comprises a system of sufficient energy to produce pions by the photoproduc- 
tion reaction 

P + 73K->A+-> * + (4) 

The energy loss of the nucleon is a significant fraction of the initial energy. 
The pion photoproduction cross section is quite large above threshold due 
to resonance production (A resonance), rising quickly to 500 mb for photon 
laboratory energies of about 0.3 GeV. Subsequent decays of the neutral and 
charged pions produced in (@]) yield high-energy photons and neutrinos. 

The center of mass energy of interactions of cosmic rays with energies 
exceeding 10 18 eV colliding with microwave-background photons is sufficient 
to generate electron-positron pairs p + j^k —> P + e+ + e ~- As a conse- 
quence, the cosmic-ray particles loose energy which leads presumably to a 
reduction of the flux or a dip in the spectrum between 10 18 and 10 19 eV 
|Hill and Schramm(1985)| |Berezinsky(2005)| |Berezinsky et al.(2004)| . 



The effect of both processes on the observed energy spectrum is frequently 
expressed by a modification factor f(E) = I p (E)/Iq{E), describing the ratio 
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Figure 2: Left: Modification factor H E) — I P (E) / I (E) of the cosmic-ray 
energy spectrum Aloisio et al.(2007a)| . Right: Loss length for protons for 
pair production and pion photoproduction |De Marco and Stanev(2 005) |. 



of the observed spectrum I p and the initial spectrum Jo as function of energy. 
The modification factor according to recent calculations is shown in Fig. [2] 
(left) [Aloisio et al.(2007a)| . A twofold structure can be recognized. A de- 
pression (the dip) at energies between 10 18 and 10 19 eV and the GZK feature 
at energies exceeding 5 • 10 19 eV. The two cases (1 and 2) represent initial 
spectra with a spectral index of 2.0 and 2.7, respectively. 

The energy loss length for pair production and pion photoproduction is 
depicted in Fig. [2] (right) |De Marco and Stanev(2005)| . Particles with ener- 
gies above the threshold of the GZK effect can travel less than about 100 Mpc 
through the Universe, before their energy has decreased to 1/e of their initial 
value. Or, in other words, particles reaching the Earth at these energies have 
propagated less than 100 Mpc, see also [Aharonian and Cronin(199 4)], and 
their sources are inside a sphere with this radius 



2.3 Multi Messenger Approach 

It has been discussed that for both scenarios, acceleration and top-down 
models, hadronic cosmic rays are accompanied by high-energy photons and 
neutrinos. Also during the propagation of hadronic particles through the 
Universe high-energy photons and neutrinos are produced. To clarify the 
origin of the highest-energy particles in the Universe, simultaneous observa- 
tions are desired of high-energy charged particles, photons, and neutrinos — 
a multi-messenger approach. Thus, the observation of high-energy charged 
particles, or charged particle astronomy, is complementary to observations 
in gamma ray astronomy [Ong(1998) Qng(2005 )] and neutr ino astronomy 



Spiering(200"3)| |Halzen(2005)[ |Berezinsky(2006)[ |Lipari(20 06) | 



Attention has to be paid on the 'simultaneous' observation: if a charged 
particle is deflected by an angle 9 in a (for this estimate simply homogeneous) 
magnetic field, its path L c h is somewhat longer as compared to the path of 
a massless neutral particle L 7 . Their relative difference can be approximated 

as 

Lch 27T6 

L 7 360°sin(e)' {) 

If a charged particle from a source at a distance L 1 = 100 Mpc is deflected by 
6 = 3°, a value 1Z = 1 + 4.6- is obtained. This corresponds to a difference 
in the arrival time of a charged particle relative to a photon (both traveling 
at the speed of light) of about AT = 150 • 10 a. Thus, for a simultaneous 
detection the acceleration processes have to be stable over such a period in 
time. 



3 Detection Method 

The extremely steeply falling cosmic-ray energy spectrum (oc E~ 3 ) yields very 
low fluxes for the highest-energy particles. At the highest energies less than 
one particle is expected per square kilometer and century. This necessitates 
huge detection areas and large measuring times. At present, they are only 
realized in huge ground based installations, registering secondary particles 
produced in the atmosphere. 

3.1 Extensive air showers 

When high-energy cosmic-ray particles penetrate the Earths atmosphere they 
interact and generate a cascade of secondary particles, the extensive air show- 
ers. Hadronic particles interact and produce new hadronic particles or gen- 
erate muons and photons through pion decays. Some of the muons may 
decay into electrons, while the photons and electrons/positrons regenerate 
themselves in an electromagnetic cascade. The by far dominant particles in a 
shower are electromagnetic particles (photons, electrons, and positrons). Most 
of the energy of the primary particle is absorbed in the atmosphere. However, 
a small fraction of the energy is transported to ground level and may be regis- 
tered in detectors for electrons, muons, and hadrons. Particles traveling with 
relativistic speeds through the atmosphere (mostly electrons and positrons) 
emit Cerenkov light. The shower particles also excite nitrogen molecules in 
the air which in turn emit fluorescence light. While the Cerenkov light is 
collimated in the forward direction of the particle, the fluorescence light is 
emitted isotropically, thus, a shower can be "viewed from aside". 

The objective of experiments observing extensive air showers is to de- 
termine the properties of the primary particle (energy Eq, mass A, arrival 
direction). In the energy regime of interest (E > 10 17 eV) mainly two 
methods are applied. Electrons (and positrons) as well as muons reach- 
ing ground level are observed in large arrays of detectors and the fluo- 
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rescence light is viewed by imaging telescopes. An alternative technique, 
presently under investigation, is the detection of radio emission from air show- 
ers. Electrons and positrons are deflected in the Earths magnetic field and 
emit synchrotron radiation, which is detected in arrays of dipole antennae 
|Falcke et al.(2005)[ [van den Berg et al.(2007)[|Huege and Falcke(2005)] . 



3.2 Measuring technique 

The direction of air showers is inferred applying two techniques. The parti- 
cles in a shower travel with nearly the speed of light through the atmosphere 
in a thin disc with a thickness of a few meters only. With detectors measur- 
ing the arrival time of the particles with a resolution of a few ns the angle 
of the shower front relative to the ground can be inferred, with the arrival 
direction being perpendicular to the shower plane. With imaging fluorescence 
telescopes the shower-detector plane is determined from the observed track in 
the camera. The orientation of the shower axis in this plane is then obtained 
by measurements of the arrival time of the photons at the detector. Using 
two (or more) telescopes to view the same shower allows a three-dimensional 
reconstruction of the shower axis. 

The shower energy is proportional to the number of electrons N e and 
muons in the shower. A simple numerical model |Horan dcl(2007bj] yields 
the relations 

E = 3.01 GeV • A om ■ N^ 96 and E = 20 GeV • A~ aA1 ■ N^ 11 (6) 

to estimate the primary energy. This illustrates that measuring N e or gives 
a good estimate for the energy almost independent of the particles mass. 
With imaging fluorescence telescopes the amount of fluorescence light can be 
measured as function of depth in the atmosphere. The total amount of light 
collected is proportional to the shower energy. Using the number of electrons 
at shower maximum, the number of photons registered per square meter in a 
detector at a distance r to the maximum of a shower with energy Eq can be 
estimated as 

7 ~ W ~ 79 V A UevJ (,/10km) 2 ' (?) 

where iV 7 w 4 7/m is the fluorescence yield of electrons in air and Xq = 
36.7 g/cm 2 (or 304 m at normal pressure) the radiation length. Absorption 
and scattering in the atmosphere have been neglected in this simple estimate, 
thus, the equation gives an upper limit for the registered photons. 

Experimentally most challenging is the estimation of the mass of the 
primary particle. Showers induced by light and heavy particles develop dif- 
ferently in the atmosphere. The depth in the atmosphere X max at which 
the showers contain a maximum number of particles depends on the primary 
particles mass 

X max = X max ~ X m A (8) 
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where is the depth of the shower maximum for proton-induced showers 

|Matthews(2005) Horandel(2007b)| . Experiments measuring the longitudinal 



shower profile by observations of fluorescence light estimate the mass by mea- 
surements of X max . 

If a shower develops higher in the atmosphere more particles (mostly elec- 
trons) are absorbed on the way to the ground. On the other hand, at high 
altitudes (with low air densities) charged pions are more likely to decay, thus, 
yielding more muons. Hence, the electron-to-muon ratio observed at ground 
level depends on the mass of the primary particle. A Heitler model of hadronic 
showers |H6randel(2007b)| yields the relation 

This implies that the registered electron-to-muon ratio depends on the energy 
per nucleon of the primary particle. 

3.3 Cosmic-Ray Detectors 

In the following we describe the most important recent detectors for ultra 
high-energy cosmic rays. 



The AGASA experiment The Akeno Giant Air Shower Array (AGASA) 
was a scintillator array located in Japan (35°N, 138°E), covering an area of 
100 km 2 [Ta keda et al.(2003)| . It consisted of 111 scintillation counters to 
register the electromagnetic shower component. Each station covered 2.2 m 2 
in area. The scintillator blocks with a thickness of 5 cm were viewed by 
a 125 mm diameter photomultiplier tube. To register the muonic shower 
component, proportional counters were used with a cross section of 10 x 10 cm 2 
and a length of 2 m or 5 m. The absorber consisted either of a 1 m thick 
concrete block, a 30 cm thick iron plate, or a 5 cm lead plate above a 20 cm 
thick iron plate. The threshold energy for muons is about 0.5 GeV. In total, 
27 detector stations were installed with areas varying from 2.8 m 2 to 20 m 2 . 



The HiRes experiment The High Resolution Fly's Eye experiment 
(HiRes) was located in Utah, USA (40°N, 112°W) |Abu-Zayyad et al.(2000a")| . 
It was the successor of the Fly's Eye experiment |Baltrusaitis et al.(19 88)|, 
which pioneered the detection of fluorescence light from air showers. HiRes 
consisted of two detector sites (Hires I & II) separated by 12.6 km, providing 
almost 360° azimuthal coverage, each. Both telescopes were formed by an 
array of detector units. The mirrors consisted of four segments and formed 
a 5.1 m 2 spherical mirror. At its focal plane an array of 16 x 16 photo- 
multiplier tubes was situated, viewing a solid angle of 16° x 16°. HiRes I 
consisted of 22 detectors, arranged in a single ring, overlooking between 3° 
and 17° in elevation. This detector used an integrating ADC read-out system, 
which recorded the photomultiplier tubes' pulse height and time information. 
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Figure 3: Schematic view of a water Cerenkov detector (left) 
and a fluorescence telescope (right) of the Pierre Auger Observatory 
|Abraham et al.(2004)| . 

HiRes II comprises 42 detectors, set up in two rings, looking between 3° and 
31° in elevation. It was equipped with a 10 MHz flash ADC system, recording 
pulse height and timing information from its phototubes. 

The Pierre Auger Observatory The observatory combines the obser- 
vation of fluorescence light with imaging telescopes and the measurement of 
particles reaching ground level in a "hybrid approach" |Abraham et al.(2004)| . 
The southern site (near Malargue, Argentina, 35.2° S, 69.5° W, 1400 m above 
sea level) of the worlds largest air shower detector is almost completed. It will 
consist of 1600 polyethylene tanks set up in an area covering 3000 km 2 . Each 
water Cerenkov detector has 3.6 m diameter and is 1.55 m high, enclosing a 
Tyvak liner filled with 12 m of high purity water, see Fig.[3j The water is 
viewed by three PMTs (9 in diameter) . Signals from the PMTs are read by the 
electronics mounted locally at each station. Power is provided by batteries, 
connected to solar panels, and time synchronization relies on a GPS receiver. 
A radio system is used to provide communication between each station and a 
central data acquisition system. 

Four telescope systems overlook the surface array. A single telescope system 
comprises six telescopes, overlooking separate volumes of air. A schematic 
cross-sectional view of one telescope is shown in Fig. [31 Light enters the bay 
through an UV transmitting filter. A circular diaphragm (2.2 m diameter), 
positioned at the center of curvature of a spherical mirror, defines the aperture 
of the Schmidt optical system. A 3.5 m x 3.5 m spherical mirror focuses the 
light onto a camera with an array of 22 x 20 hexagonal pixels. Each pixel has 
a photomultiplier tube, complemented by light collectors. Each camera pixel 
has a field of view of approximately 1.5°, a camera overlooks a total field of 
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Figure 4: Accumulated ex- 
posures of various experi- 
ments |Kampert(20Q8)| . 

view of 30° azimuth x 28.6° elevation. 

Telescope Array Like the Pierre Auger Observatory, the Telescope array 
is a hybrid detector, presently under construction in Mullard County, Utah, 
USA |Fukushima et al.(2007)| . It covers an area of 860 km 2 and comprises 
576 scintillator stations and three fluorescence detector sites on a triangle with 
about 35 km separation, each equipped with twelve fluorescence telescopes. 

Accumulated exposures (i.e. experiment aperture times measuring time) for 
various high-energy experiments are presented in Fig.[H [Kampert(2008)| , For 
surface arrays the aperture is a function of the detector area and constant with 
energy. On the other hand, the aperture of fluorescence detectors depends on 
the shower energy, low energy showers can be seen up to a restricted distance 
only. This may be illustrated using the approximation ([7]) : the fiducial volume 
to register z/™ m photons can be estimated as 

V fid OC (j,™™)- 1 - 5 ,4-0.069 £1.569 ( 1Q ) 

This shows that the fiducial volume is a function of the primary energy. In this 
simple approximation there is a small dependence on the mass of the primary 
particle (« 25% difference between proton and iron induced showers) and an 
increase of almost a factor of 40 in the fiducial volume per decade in primary 
energy. A similar energy dependence can be recognized in Fig.|4]for the various 
fluorescence detectors. For fluorescence telescopes with a limited field of view 
in elevation an additional effect occurs: low energy showers penetrate less 
deep into the atmosphere and thus may have their maximum above the field 
of view of the telescopes, thus, reducing further the effective aperture. 
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Figure 5: All-particle en- 
ergy spectra measured 
by the Pierre Auger Ob- 
servatory using different 
reconstruction methods 
|Yamamoto et al.(2007)]. 



The Pierre Auger Observatory, still under construction, is already the 
largest cosmic ray detector, the accumulated data exceed the data of all previ- 
ous experiments. In particular, those of the largest scintillator array (AG AS A) 
and the largest pure fluorescence detector (HiRes). Thus, the Pierre Auger 
Observatory is expected to measure the properties of ultra high-energy cosmic 
rays with unprecedented accuracy. 



4 Energy Spectrum 

Measurements of the energy spectrum provide important information about 
the origin of cosmic rays. Over a wide range in energy the all-particle differen- 
tial energy spectrum is usually described by a power law dN/dE oc E 1 . For 
energies below 10 15 eV a value for the spectral index 7 = —2.7 has been estab- 
lished by many experiments. The most prominent feature in the all-particle 
spectrum is the so called knee at an energy of about 4-10 15 eV. At this energy 
the spectral index changes to 7 w —3.1. The knee in the all-particle spectrum 
is caused by the subsequent cut-offs (or knees) of the spectra of individual 
elements, starting with protons at ss 4.5 • 10 15 eV. However, this feature 
is below the focus of the current article, thus, the reader may be referred to 
e.g. [Horandel(2003a)| |Horandel(2004)| |Horandel(2008)] for a more detailed 



discussion about galactic cosmic rays and the knee. In the following we focus 
on energies above 10 17 eV. 

Recent energy spectra as obtained by the Pierre Auger Observatory are 
depicted in Fig. [5] [Yamamoto et al.(2007)| . The registered flux has been mul- 
tiplied by E 3 . Different methods are applied to reconstruct the spectra. The 
first method uses the data from the 3000 km 2 surface array. The detection 
efficiency reaches 100% for showers with zenith angles less than 60° for ener- 
gies above 10 18 5 eV and for inclined showers (6 > 60°) above 10 18 8 eV. The 
signal at 1000 m from the shower axis is used to estimate the shower energy. 
To avoid a dependence on interaction models used in air shower simulation 
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Table 1: Energy shifts applied to individual experiments as shown in Fig.[6l 



Experiment 



Reference 



Takeda et al.(2Q03)| 



Energy shift 



AGASA 

Akeno 1 km 2 

Akeno 20 km 2 

Auger 

Fly's Eye 

Haverah Park 

HiRes-I 

HiRes-II 

HiRes-MIA 

KASCADE-Grande 

MSU 

SUGAR 

Yakutsk T500 

Yakutsk T1000 



Nagano et al.(1984a) 
Nagano et al.(1984b) 



Yamamoto et al (2007)] 
Bird et al.(1994)| 
Ave et al. (2003)1 
Abbasi et al.(2004a)| 
Abbasi et al.(2005)| 



Abu-Zayyad et al. (2 000b) | 
Haungs et al.(2008)| 
Fomin et al.(1991)] 



Anchordoqui and Goldberg(2004)| 



Glushkov et al.(2003) 



Glushkov et al.(2003) 



-22% 
-4% 
-22% 
-20% 
-3% 
-2% 
0% 
0% 
+5% 
-7% 
-5% 
0% 
-35% 
-20% 



codes, an energy estimator is derived based on measured showers: a subset 
of showers contains so called hybrid events, seen simultaneously by the sur- 
face detector array and at least one fluorescence telescope. The fluorescence 
telescopes provide a nearly model independent calorimetric energy measure- 
ment of the showers in the atmosphere. Only a small correction for 'invisible 
energy' (high-energy muons and neutrinos carrying away energy) has to be 
applied. This factor amounts to about 10% and contributes with about 4% 
to the systematic error for the energy. The energy calibration thus obtained 
is applied to all events recorded with the surface detector array. Also inclined 
events with zenith angles exceeding 60° have been analyzed, yielding the sec- 
ond spectrum displayed. Finally, a set of showers which have been recorded 
by at least one surface detector tank and one fluorescence telescope have been 
analyzed. The resulting energy spectrum reaches energies as low as 10 18 eV, 
as can be inferred from Fig.OD It is interesting to point out that the different 
spectra have been analyzed independently and agree quite good with each 
other. 

The all-particle energy spectra as obtained by various experiments are 
compiled in Fig.[6l The flux has been multiplied by E 3 . The upper panel 
shows the original data. The different experiments yield absolute values which 
differ by almost an order of magnitude in this representation. Nevertheless, 
the overall shape of the energy spectrum seems to be reflected in all data, 
irrespective of the absolute normalization. This becomes more obvious when 
the energy scales are slightly readjusted. Typical systematic uncertainties for 
the energy scale are of order of 10% to 30% in the region of interest. When 
energy shifts are applied, the results have to be treated with care since the 
apertures of some experiments change as function of energy (see Fig. [4]) and 
this effect has not been taken into account in the procedure used here. 
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Figure 6: All-particle energy spectra as obtained by different experiments. 
The top panel shows the original values, in the bottom panel the energy 
scales of the individual experiments have been adjusted. For references and 
energy shifts, see Table [U The lines indicate the end of the galactic com- 
ponent according to the poly-gonato model |Horandel(2003a)] and a possible 
contribution of extragalactic cosmic rays. 
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In the lower panel of Fig. [6] the energy scales of the different experi- 
ments have been adjusted to fit the flux according to the poly-gonato model 
at 10 18 eV. The latter has been obtained through a careful procedure ex- 
trapolating the measured spectra for individual elements at low energies 
|Horandel(2003a)| . Thus, the normalization applied provides a consistent de- 
scription from direct measurements (10 GeV region) up to the highest energies. 
The corresponding energy shifts are listed in Table [TJ 

The normalized spectra agree very well and seem to exhibit a clear shape 
of the all-particle energy spectrum. Some structures seem to be present in 
the spectrum. The second knee at about 4 • 10 17 eV, where the spectrum 
steepens to 7 ~ —3.3 and the ankle at about 4 • 10 18 eV, above this energy 
the spectrum flattens again to 7 ~ —2.7. Finally, above 4 • 10 19 eV the 
spectrum exhibits again a steepening with a spectral index 7 w —4 to —5. 
The new Auger results help to clarify the situation in this energy region. 
While the AGASA experiment has reported events beyond the GZK threshold 
|Take da et al.(1998)], the HiRes e xperiment has reported a detection of the 
GZK cut-off [Abbasi et al.(2007)| . With the new results, a steeper falling 
spectrum above 4 • 10 ly eV is now confirmed. 

The second knee possibly marks the end of the galactic component 
[Horandel(20Q3a)| . If the energy spectra for individual elements exhibit knees 
at energies proportional to their nuclear charge, the heaviest elements in galac- 
tic cosmic rays should fall off at an energy of about 92 • E% w 4 • 10 17 eV. An 
interesting coincidence with the position of the second knee. Different scenar- 
ios for the transition from galactic to extragalactic cosmic rays are discussed 
e.g. in |Kampert(2007) H6randel(2008)|. In the energy region around the 
ankle a depression is seen in the all-particle flux, also referred to as the dip. 
It is proposed that this dip is caused by interactions of ultra high-energy par- 
ticles with the cosmic microwave background, resulting in electron-positron 
pair production, see Sect. 12.21 The steepening in the flux above 4 • 10 19 eV 
could be an indication of the GZK effect, i.e. photo-pion production of ultra 
high-energy cosmic rays with the microwave background, see Sect. 12.21 How- 
ever, for a definite answer also other properties of cosmic rays have to be 
investigated. 



5 Mass Composition 

The elemental composition of galactic cosmic rays has been discussed else- 
where, e.g. |Horandel(2003b)[|Horandel(2008)| . Above 10 17 eV the situation 



is experimentally very challenging, since we are far away in parameter space 
from collider experiments, where the properties of high-energy interactions 
are studied in detail. Thus, the air shower models used to interpret the data 
have to extrapolate over a wide range in parameter space. 

The fraction of iron nuclei in cosmic rays as deduced by many experiments 
has been investigated |Dova et al.(2005)] . No clear conclusion can be drawn 
about the composition at highest energies. Tension in the interpretation of 
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the measured data has been observed as well by the HiRes-MIA experiment 
|Abu-Zayyad et al. (2000c)] . The observed X. max values exhibit a trend to- 
wards a lighter composition as function of energy in the range between 10 17 
and 10 18 eV. On the other hand, measured muon densities indicate a very 
heavy composition in the same energy range. 

Methods relying on the measured muon densities, the lateral distribution 
of Cerenkov light registered at ground level, or geometrically-based methods 
are rather indirect and depend on certain assumptions and/or interaction 
models. The most bias free mass estimator is probably a measurement of the 
depth of the shower maximum X max , preferably with an imaging telescope 
such as fluorescence detectors. The best way to infer the mass is to measure 
X max distributions, rather than average values only. However, unfortunately, 
also the interpretation of the measured values depends on hadronic interaction 
models used in air shower simulations. 

The average depth of the shower maximum registered by several experi- 
ments is plotted in Fig.[7|as function of energy. In the top panel the data are 
compared to predictions of air shower simulations for primary protons and 
iron nuclei, using different hadronic interaction models, namely QGSJET 01 
Kalmykov et al.(1997)] , QGSJET II-3 |Ostapchenko(2005)1 , SYBILL 2.1 
Engel et al.(1999)| , and DPMJET 2.55 |Ranft(1995)| . The models yield dif- 
ferences in X max of order of 30 g/cm 2 for iron nuclei and w 50 g/cm 2 for 
proton induced showers. An overall trend seems to be visible in the data, 
the measured values seem to increase faster with energy as compared to the 
model predictions. This implies that the composition becomes lighter as func- 
tion of energy. Through interactions with the cosmic microwave background 
heavy nuclei are expected to break up during their propagation through the 
Universe (GZK effect) and a light composition is expected at the highest en- 
ergies. However, e.g. the Auger data at the highest energies correspond to a 
mixed composition for all models displayed. 

In the bottom panel of Fig. [7] the measured values are compared to pre- 
dictions of astrophysical models of the origin of high-energy cosmic rays. 
The propagation of high-energy cosmic rays in extragalactic turbulent mag- 
netic fields is considered in |Globus et al.(2007)| . The average X max values 
are shown for a case, assuming a mixed source composition with an injection 
spectrum cx E~ 2A , a continuous distribution of the sources, and no extra- 
galactic magnetic field. Other cases studied deliver similar results in X max , 
for details see |Globus et al.(2007)| . 

Different scenarios for the transition from galactic to extragalactic cosmic 
rays are discussed in |Aloisio et al.(2007b)| . Two scenarios are distinguished, 
a 'dip' model in which the galactic and extragalactic fluxes equal at an en- 
ergy below 10 18 eV and an 'ankle' approach in which both components have 
equal fluxes at an energy exceeding 10 19 eV. It is proposed that the dip is a 
consequence of electron-positron pair production, see Sect. 12.21 For the 'dip' 
model a source spectrum cx E~ 2 ' 7 is assumed and a spectrum cx E~ 2 for the 
ankle approach. Cosmic rays have been propagated through an extragalactic 
magnetic field of 1 nG. The resulting average X max values, based on simula- 
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Figure 7: Average depth of the shower maximum X max as function of energy 
as measured by the Pierre Auger Observatory [Unger et al.(2007)] , as well as 
the Fly's Eye [Bird et al.(1994)| , H averah Park [Watson(2000)| , HiRes/MIA 
Abu-Zayyad et al.(2000c)| , HiRes |Abu-Zayyad et al.(2000b)| , and Yakutsk 
Knurenko et al.(2001)| experiments. Top: measured values are compared 
to predictions for primary protons and iron nuclei for different hadronic in- 
teraction models QG SJET 01 |Kalmykov et al.(1997)| (— ), QGSJET II-3 
|Ostapchenko(20Q5)| (- --), SYBILL 2.1 [Engel et al.(1999)| (•••), and DPM- 
JET 2.55 [Ranft(1995)| (■ — ■). B ottom: comparison to astrophysical models 
according to |Globus et al.(2007j] ( — ) as well as |Aloisio et al. (2007b j] , for 
the latter a dip ( — ) and an ankle (• • •) scenario are distinguished. 



16 




EJeV] 



Figure 8: Upper limits on the fraction of photons in the integral cosmic- 
ray flux compared to predictions for GZK photons and top-down scenar- 
Abraham et al.(2007a)|. Experimental da ta are from the Auger sur- 



face deTeT-torTTarrows^ and a hybrid analysis (FD) 

|Abraham et al.(2007b)|, Haverah Park (HP ) [Ave et al. (2002)1 , AG AS A 
(A) |Shinozaki et al. (2002)1 |Risse et al.(20Q5)] AG AS A and Yakuts k (AY) 
|Rubtsov et al. (2006)| , as well as Yakutsk (Y) |Glushkov et al.(2007)] . 



tions using the interaction code QGSJET 01 are displayed in the figure. 
The figure illustrates that we are entering an era where it should be possible 
to distinguish between different astrophysical scenarios. 

Of great interest is also whether other species than atomic nuclei contribute 
to the ultra high-energy particle flux. 



5.1 Photon Flux Limit 

Air showers induced by primary photons develop an almost pure electro- 
magnetic cascade. Experimentally they are identified by their relatively low 
muon content or their relatively deep shower maximum. Since mostly elec- 
tromagnetic processes are involved in the shower development, the predic- 
tions are more reliable and don't suffer from uncertainties in hadronic in- 
teraction models. A compilation of recent upper limits on the contribution 
of photons to the all-particle flux is shown in Fig. [8] [Abraham et al.(2007a)| . 
The best photon limits are the latest results of the Pierre Auger Observa- 
tory Abraham et al.(2007a)| setting rather strong limits on the photon flux. 
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Figure 9: Limits at 90% confidence level for a diffuse flux of v T assuming 
a 1:1:1 ratio of the three neutrino flavors at Earth |Abraham et al.(2007c)| 
Kampert (2008)| . The experimental results are compared to predictions for 



GZK neutrinos and a top-down model Kalashev et al.(2002)|. 



They are based on measurements with the Auger surface detectors, taking 
into account observables sensitive to the longitudinal shower development, 
the signal rise time, and the curvature of the shower front. The photon frac- 
tion is smaller than 2%, 5.1%, and 31% above energies of 10 19 , 2 • 10 19 , and 
4 • 10 19 eV, respectively with 95% confidence level. 

In top-down scenarios for high-energy cosmic rays, the particles are decay 
products of super-heavy particles. This yields relatively high-fluxes of pho- 
tons predicted by such models. Several predictions are shown in the figure 
|Aloisio et al.(2004)[ |Ellis et al.(2006)| . These scenarios are strongly disfa- 



vored by the recent Auger results. 

The upper limits are already relatively close to the fluxes expected for 
photons originating from the GZK effect |Gelmini et al.(2005)| , shown in the 
figure as shaded area. 

5.2 Neutrino Flux Limit 

The detection of ultra high-energy cosmic neutrinos is a long standing ex- 
perimental challenge. Many experiments are searching for such neutrinos, 
and there are several ongoing efforts to construct dedicated experiments to 
detect them |Halzen et al.(2002)[ [Falcke et al.(2004)| . Their discovery would 



open a new window to the Universe Becker(2007)| . However, so far no ultra 



high-energy neutrinos have been detected. 



2 



2 Neutrino s produced in air showers (atmospheric neutrinos) Fu kuda et al.(1998)_ 
in the sun |Fukuda et al.(1999)[ [Ahmed et al.(2004)| , and during super nova 1987A 
[Hirata et al.(1987)[|Bionta et al.(1987)| have been detected, but are at energies much below 
our focus. 
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As discussed above (Sect. [2]) ultra high-energy cosmic rays are expected to 
be accompanied by ultra high-energy neutrinos. The neutrinos are produced 
with different abundances for the individual flavors, e.g. pion decay leads to a 
ratio v e : = 2 : 1. However, due to neutrino oscillations the ratio expected 
at Earth is v T : : u e = 1 : 1 : 1. 

To discriminate against the huge hadronic background in air shower de- 
tectors, neutrino candidates are identified as nearly horizontal showers with 
a significant electromagnetic component. The Pierre Auger Observatory is 
sensitive to Earth-skimming tau-neutrinos that interact in the Earth's crust. 
Tau leptons from v T charged-current interactions can emerge and decay in 
the atmosphere to produce a nearly horizontal shower with a significant elec- 
tromagnetic component. Recent results from the Pierre Auger Observatory 
together with upper limits from other experiments are presented in Fig.[9l 
Assuming an E~ 2 differential energy spectrum Auger derives a limit at 90% 
confidence level of E„ AN Vt /6E v < 1.3 -10~ 7 GeV cm -2 s _1 sr" 1 in the energy 
range between 2 • 10 17 and 2 • 10 19 eV. 

According to top-down models for ultra high-energy cosmic rays a large 
flux of ultra high-energy neutrinos is expected. As an example, the predic- 
tions of a model |Kalashev et al.(2002)| are shown in the figure as well. This 
model is disfavored by the recent upper limits. It should also be noted that 
the current experiments are only about one order of magnitude away from 
predicted fluxes of GZK neutrinos (cosmogenic neutrinos). 



6 Arrival Direction 

The arrival directions of cosmic rays provide an important observable to 
investigate the sources of these particles. Since charged particles are de- 
flected in magnetic fields, the cosmic-ray flux observed at Earth is highly 
isotropic. A significant evidence for an anisotropy in the arrival directions 
would be the most direct hint towards possible cosmic-ray sources. Unfortu- 
nately, only limited experimental information is available about both, galactic 
Rand and Kulkarni(1989)[ |Vallee(2004)| and extragalactic [Kronberg(1994) 



Grasso and Rubinstein(200T)| magnetic fields. Selecting particles at the high 



est energies limits the field of view to distances less than 100 Mpc, see Fig. [2] 
(right). This implies two advantages: the number of source candidates is 
limited and the particles are only slightly deflected since they propagate a 
restricted distance only. 



6.1 Galactic Center 

The center of our galaxy is an interesting target for cosmic-ray anisotropy 
studies. It harbors a massive black hole associated with the radio source 
Sagittarius A* and a supernova remnant Sagittarius A East. Both are candi- 
dates to be powerful cosmic-ray accelerators. The importance is underlined 
by recent discoveries: the HESS experiment has reported the observation from 
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Figure 10: Significance maps of excess/deficit events in equatorial coordi- 
nates as measured by the AGASA Hayashida et al.(1999)| (left) and SUGAR 
Belli do et al.(2001)] (right) experiments. The lines in both panels indicate 
the galactic plane. AGASA: events within a radius of 20° are summed up in 
each bin. SUGAR: The white circle with a radius of 5.5° indicates the error 
for a point source. 




Figure 11: Map of cosmic-ray 
over-density significances 
near the Galactic Cen- 
ter as measured by the 
Pierre Auger Observatory 
|Abraham et al.(2007d)| . The 
line indicates the galactic 
plane and " + " marks the 
Galactic Center. The circles 
represent the regions of excess 
events seen by the AGASA 
and SUGAR experiments. 



TeV 7 rays near the location of Sagittarius A* |Aharonian et al.(2004)| and 
discovered a region of extended emission from giant molecular clouds in the 
central 200 pc of the Milky Way [Aharonian et al.(2006)] . 

Of particular interest to search for anisotropies in cosmic rays is the region 
around 10 18 eV. At these energies the tail of the galactic component might 
still contribute significantly to the all-particle spectrum and neutrons from 
the galactic center can reach the Earth without decaying. Such neutrons 



would not be deflecte d by magnetic fields [Medina Tanco and Watson(2001) 
Bossa et al.(2003)[ |Aharonian and Neronov(2005)[ |Crocker et al.(2005) 



Grasso and Maccione(2005) Biermann et al.(2004) 

The AGASA experiment has investigated anisotropies in the ar- 
rival directions of cosmic rays at energies around 10 18 eV, see Fig.fTOl 
|Hayashi da et al.(1999)| . The Galactic Center is just outside the field of view 
of the experiment. However, an excess in the Galactic-Center region has been 
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Figure 12: Arrival directions of cosmic rays with energies exceeding 4- 10 19 eV 
in equatorial coordinates as observed by the AG AS A experiment. Red squares 
and green circles represent cosmic rays with energies exceeding 10 20 eV and 
(4 - 10) • 10 19 eV, respectively. |Takeda et al.(1999)||AGASA(2003)| 



detected. Also the SUGAR experiment, located in Australia has reported an 
excess of events from the region of the Galactic Center at 10 18 eV, see Fig.fTOl 
(right) [Bellido et al.(2001)| . It should be noted that both findings are on the 
3 to 4cr level only. 

Recently, data from the Pierre Auger Observatory have been searched for 
anisotropies in the region of the Galactic Center |Abraham et al.(2007d)| . A 
map of resulting cosmic-ray over-density significances is displayed in Fig.fTTl 
The regions were AGASA and SUGAR have found an excess are marked in the 
figure. With a statistics much greater than those of previous experiments, it 
has been searched for a point-like source in the direction of Sagittarius A. No 
significant excess has been found. Also searches on larger angular scales show 
no abnormally over-dense regions. These findings exclude recently proposed 
scenarios for a neutron source in the Galactic Center. 



6.2 Clustering of Arrival Directions 

The AGASA experiment has investigated small-scale anisotropies in the ar- 
rival directions of cosmic rays Takeda et al.(1999)] . Above an energy of 



4 • 10 19 eV they have found clusters of events coming from the same direc- 
tion, see Fig. [321 One triplet and three doublets with a separation angle of 
2.5° have been reported, the probability to observe these clusters by a chance 
coincidence under an isotropic distribution is smaller than 1%. 

The HiRes experiment has found no significant clustering at any angular 
scale up to 5° for energies exceeding 10 EeV Abbasi et al. (2004b)] . Combining 



data from the AGASA, HiRes, SUGAR, and Yakutsk experiments at energies 
above 40 EeV a hint for a correlation has been found at angular scales around 
25° |Kachelriess and Semikoz(2006)| . 

Also the data of the Pierre Auger Observatory have been searched for 
clustering in the arrival directions |Mollerarch et al.(2007)| . The autocor- 
relation function has been analyzed adopting a method, in which a scan 
over the minimum energy E and the separation angle O is performed 
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Figure 13: Autocorre- 
lation scan for events 
with energies above 
10 20 eV recorded with the 
Pierre Auger Observatory 
|Mollerarch et al.(20Q7j] . 
The chance probability is 
shown as function of sepa- 
ration angle and threshold 
energy. 



[Finley and Westerhoff(2004)| . For each value of E and a chance prob- 
ability is calculated by generating a large number of isotropic Monte Carlo 
simulations of the same number of events, and computing the fraction of sim- 
ulations having an equal or larger number of pairs than the data for those 
parameters. The result is depicted in Fig.rjJl showing the probability as func- 
tion of separation angle and threshold energy. A broad region with an excess 
of correlation appears at intermediate angular scales and large energies. The 
minimum is found at 7° for the 19 highest events (E > 57.5 EeV), where 
eight pairs were observed, while one was expected. The fraction of isotropic 
simulations with larger number of pairs at that angular scale and for that 
number of events is P m in — 10 -4 . The chance probability for this value to 
arise from an isotropic distribution is P w 2 • 10 -2 . 



6.3 Correlation with BL-Lacs 

Interesting candidates as cosmic-ray sources are BL Lacertae objects. They 
are a sub class of blazars, active galaxies with beamed emission from a rela- 
tivistic jet which is aligned roughly towards our line of sight. Several exper- 
iments have searched for correlations of the arrival directions of cosmic rays 
with the position of BL Lacs on the sky. 

A correlation was found between a subset of BL Lac positions 
and arrival directions recored by AGASA with energies exceeding 
48 EeV and by the Yakutsk experiment at energies above 24 EeV 
|Tinyakov and Tkachev(2001)|. This correlation and further ones as re- 



ported in |Tinyakov and Tkachev(2002)[ |Gorbunov et al.(2002)| between BL 
Lacs and ultra high-energy cosmic rays registered by the AGASA and 
Yakutsk experiments were not confirmed by data of the HiRes experiment 
|Abbasi et al.(2006)| . On the other hand, an excess of correlations was 
found for a subset of BL Lacs and cosmic rays with energies above 10 EeV 
|Abbasi et al.(2006)[|Gorbunov et al.(2004)| . 

In spring 2007 the number of events recorded by the Pierre Auger Ob- 
servatory above 10 EeV was six times larger than the data used in previous 
searches. The correlation hypotheses reported previously have been tested 
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Figure 14: Number of events correlated with confirmed BL Lacs with optical 
magnitude m < 18 from the 10 th edition of the catalog of quasars and nu- 
clei |Veron-Cetty and Veron(2006)| (points) and average for an isotropic flux 
(solid line) along with dispersion in 95% of simulated isotropic sets (bars) 
|Har ari et al.(2Q07)] . As function of the angular separation (for E > 10 EeV, 
left) and as function of threshold energy (0 < 0.9°, right). 



with the Auger data Harari et al.(2007j] . Since the southern detector of the 
Pierre Auger Observatory sees a different part of the sky as compared to 
the AGASA, HiRes, and Yakutsk experiments, only the 'recipes' could be 
tested but using different sources on the sky. Non of the previously reported 
hypotheses could be confirmed, the chance probabilities for the different ap- 
proaches were found to be slightly smaller than 1%. The correlations search 
has been extended to a broader range of angular scales and energy thresholds, 
see Fig.HH It shows the number of correlated events as function of separation 
angle (left) and energy threshold (right). The curves represent expectations 
for an isotropic flux. The error bars depict the dispersion within 95% of sim- 
ulated isotropic sets. As can be inferred from the figure, the measured data 
are compatible with an isotropic distribution and they do not confirm earlier 
findings. 



6.4 Correlation with AGN 

Another interesting set of objects to serve as sources of ultra high-energy cos- 
mic rays are Active Galactic Nuclei (AGN). The radiation from AGN is be- 
lieved to be a result of accretion on to the super-massive black hole (with 10 6 to 
10 s solar masses) at the center of the host galaxy. AGN are the most luminous 
persistent sources of electromagnetic radiation in the Universe. An example 
of an AGN is shown in Fig.|T5j Centaurus A is with a distance of 3.4 Mpc one 
of the closest AGN. The radio lobes are thought to be the result of relativistic 
jets emerging from the central black hole. Different scenarios related to AGN 
have been developed, which are supposed to accelerate particles to highest en- 
ergies, e.g. [Biermann and Strittmatter(1987) Rachen and Biermann(1993) , 



Romero et al.(1996)[ |Ostrowski(1999)[ |Lyutikov and Ouyed(1 999)| 



The arrival directions of cosmic rays as measured by the Pierre 
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Figure 15: Centaurus A 
as seen by the Hub- 
ble Space Telescope and 
the VLA radio telescope 
|http://hubblesite.org| . The 
radio lobes extend over a 
scale of about 10° along the 
super-galactic plane. 




Figure 16: Aitoff projection of the celestial sphere in galactic coordinates 
with circles of 3.2° centered at the arrival directions of 27 cosmic rays 
detected by the Pierre Auger Observatory with energies E > 57 EeV 
|Abraham et al.(2007e)] . The positions of AGN with redshift z < 0.018 
(D < 75 Mpc) from the 12 th edition of the catalog of quasars and nuclei 
|Veron-Cetty and Vero n(2006)| are indicated by the asterisks. The solid line 
draws the border of the field of view of the southern observatory (for zenith 
angles < 60°). Darker color indicates larger relative exposure. The dashed 
line indicates the super-galactic plane. Centaurus A, one of the closest AGN 
is marked in white. 
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Auger Observatory have been correlated with the positions of AGN 
|Abraham et al.(2007e) Abraham et al.(2007f)| . Data taken during the con- 



struction of the observatory since January 2004 have been analyzed, corre- 
sponding to slightly more than one year of data of the completed observatory. 
The angular resolution of the detector is better than 1° at energies above 
10 19 eV |Ave et al.(2007a)| . The positions of AGN according to the 12 th 
edition of the catalog of quasars and nuclei |Veron-Cetty and Veron(2 006) 
within a distance D have been used. A scan has been performed over the 
distance D, a threshold energy E t h, and the correlation angle O. The best 
correlation has been found for events with energies exceeding E t h — 57 EeV, 
a maximum distance D = 71 Mpc, corresponding to a redshift z = 0.017, and 
a correlation angle = 3.2°. With these parameters 20 out of 27 cosmic rays 
correlate with at least one of the 442 selected AGN (292 in the field of view of 
the observatory). Only 5.6 are expected, assuming an isotropic flux. The 27 
cosmic rays measured with the highest energies are shown in Fig. [TBI together 
with the positions of the AGN. Many of the observed correlated events are 
aligned with the super-galactic plane. Two events have arrival directions less 
than 3° away from Centaurus A. These results indicate clearly that the arrival 
directions of cosmic rays at highest energies are not isotropic. 

A cosmic ray with charge Ze that travels a distance D in a regular magnetic 
field B is deflected by an angle Abraham et al.(2007f)| 



60EeV 

E/Z 



D ( dx B 



x 



\kpc 3 /iG 



(11) 



Assuming a coherence scale of order « 1 kpc |Stanev(1997)| for the regular 
component of the galactic magnetic field, the deflection angle is a few degrees 
only for protons with energies larger than 60 EeV. This illustrates that the 
observed angular correlations are reasonable, but one has to keep in mind 
the limited knowledge about galactic magnetic fields. The angular scale of 
the observed correlations also implies that intergalactic magnetic fields along 
the line of sight to the sources do not deflect cosmic-ray trajectories by much 
more than a few degrees. The root mean square deflection of cosmic rays with 
charge Ze, traveling a distance D in a turbulent magnetic field with coherence 
length L c is [Abraham et al.(2007f)| 



_ AO &Q EeVB rms I D I L c 
d ™ s E/Z nG Vl00Mpc\/lMpc- [ ' 

As information on intergalactic magnetic fields is very sparse, the correlations 
observed can be used to constrain models of turbulent intergalactic magnetic 
fields. Within the observed volume they should be such that in most directions 
B rms \fL~c < 10~ 9 G\/Mpc. In the future the Pierre Auger Observatory will 
collect more data and more than one event per source should be detected. It 
should then be possible to use the data itself to set constraints on magnetic 
field models. 
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It should be noted that the findings by the Pierre Auger Observatory imply 
that the sources of the highest-energy-cosmic rays are spatially distributed like 
AGN. The actual acceleration sites could be the AGN itself or other candidates 
with the same spatial distribution as AGN. 



7 Discussion and Outlook 

"How do cosmic accelerators work and what are they accelerating?" is one 
of eleven science questions for the new century in physics and astronomy 
[National Research Council(2003)| . In the last few years important progress 
has been made in measuring the properties of ultra high-energy cosmic rays. 
In particular, the results of the Pierre Auger Observatory have significantly 
contributed to an improvement in understanding the origin of the highest- 
energy particles in the Universe. The discovery of correlations between the 
arrival directions of cosmic rays and the positions of AGN was among the 
most important scientific breakthroughs in 2007 for several science media 
organizations, see [www.auger.org/news/top_news_2007.html]. 

The must important findings discussed in this overview may be summa- 
rized as follows. The structures in the energy spectrum at highest energies 
seem to become more clear. In particular, there seems to be evidence for a 
steeper falling spectrum above 4 • 10 19 eV (Figs.[5]and_]). The question arises 
whether this steepening is due to the GZK effect or due to the maximum 
energy achieved during the acceleration processes. The most convincing evi- 
dence for the existence of the GZK effect is provided by the correlations of the 
arrival directions with AGN. They occur sharply above an energy of 57 EeV. 
At this energy, the flux measured by the Pierre Auger Observatory is about 
50% lower than expected from a power law extrapolation from lower energies, 
see Fig.[5l Thus, there seems to be a connection between the steepening in 
the spectrum and the AGN correlation. 

The correlations occur on an angular scale of about 3.2°. This indicates 
that the particles are deflected marginally only. In turn, this implies they 
should be light particles, with low Z , see ifTTj) and (fl2|) . However, there is some 
tension between this expectation and the measurements of the average depth 
of the shower maximum X max (Fig.[7j). The data at 4- 10 19 eV are compatible 
with a mixed composition. But, since the correlations occur relatively sharp 
above 57 EeV, some dramatic change in composition above this energy can 
also not be excluded. 

The correlation implies that the sources of ultra high-energy particles are 
in our cosmological neighborhood (D < 75). The GZK horizon, defined as 
the distance from Earth which contains the sources that produce 90% of the 
protons that arrive with energies above a certain threshold is 90 Mpc at 80 EeV 
and 200 Mpc at 60 EeV [Harari et al.(2006)| . There seems to be a slight 



mismatch between these numbers and the Auger findings. Shifting upward 
the Auger energy scale by about 30%, as indicated by some simulations of the 
reconstruction procedures |En gel et al.(2007)] , a better agreement between 
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the predicted GZK horizon and energy threshold with the observed data could 
be achieved Abraham et al.(2007f)|. 



The biggest uncertainty in the absolute energy scale of the Pierre Auger 
Observatory is the knowledge of the fluorescence yield. At present, intensive 
efforts are conducted by various groups to precisely determine the fluorescence 
yield of electrons in air Arqueros et al.(2Q08)| . Attention is paid to the de- 



pendence of the yield on atmospheric parameters, like pressure, temperature, 
and humidity. In particular, upcoming results from the AIRFLY experiment 
|Ave et al. (2007b) Privitera et al.(2007)J are expected to significantly reduce 



the uncertainties of the energy scale for fluorescence detectors. 

The correlation between the arrival directions and the positions of 
AGN sets constraints on models for the acceleration of ultra high-energy 
particles. The results imply that the spatial distribution of sources is 
correlated with the distribution of AGN. Thus, already some scenarios 
are strongly disfavored. Ruled out are models proposing sources in our 
Galaxy, like neu tron stars [Blasi et al.(2000)| , pulsars Bednarek(2003)|, 



and black holes [Par and~P laga(1999)|. Models for sources in the galac- 
tic halo are also ruled out such as top-down scenarios with decaying 
super-heavy particles [Berezinsky et al.(1997) Kuz'min and Rubakov(1998), 



Birkel and Sarkar(1998)| . These models are also severely constraint by the 



upper limits on the photon flux (Fig.E} and the neutrino flux (Fig. [9]). Within 
the next years measurements of photons and neutrinos produced in the GZK 
effect seem to be in reach. Their detection would be an important and com- 
plementary information about the origin and propagation of ultra high-energy 
cosmic rays. 

In summary, the acceleration of ultra high-energy particles in 
AGN seems to be very attractive, different scenarios have been pro- 



posed, e.g. [Biermann and Strittmatter(1987) Rachen and Biermann(1993) 
Romero et al.(1996)[|Ostrowski(1999)||Lyutikov and Ouyed(1999)| . However, 



other sources with a similar spatial distribution are not excluded. 

With the energy density estimated in Sect. 12.11 we obtain a total cosmic 
ray power of about 9.7 • 10 43 erg/s within a sphere (r = 75 Mpc) seen by 
the Pierre Auger Observatory at the highest energies. The typical power in 
the jets of AGN is of order of 10 44 to 10 46 erg/s |Kording et al.(2008)| . If we 
assume about 10% of this power is converted into cosmic rays, about 1 to 10 
sources are needed to sustain the power of the observed extragalactic cosmic 
ray flux within a distance of 75 Mpc from Earth. If the efficiency is slightly 
smaller, the number of sources required is correspondingly slightly larger. 3 
If the sources of the highest-energy particles are indeed related to AGN, the 
number of correlated events seen by the Pierre Auger Observatory seems to be 
of the right order of magnitude and one expects to see in future more events 
from the same sources. 

When the number of correlated events found in the Auger data is compared 
to expectations for the AG AS A and HiRes experiments, one has to be aware 



3 Based on statistical arguments a minimum number of sources > 61 has been estimated 
[Abraham et al. (20071)] . 
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of the different energy scales, see Table [TJ The energy scales of the AG AS A 
and HiRes experiment are shifted relative to the Pierre Auger Observatory 
by about 42% and 20%, respectively. If the Auger prescription is applied 
to the data of these experiments, the energy threshold has to be adjusted 
correspondingly. 

In the next years several experiments focus on the exploration of the en- 
ergy region of the transition from galactic to extragalactic cosmic rays (10 17 to 
10 18 eV). The 0.5 km 2 K ASCADE-Grande exper iment |Navarra et al.(2004)| 
is taking d ata since 2004 |Chiavassa et al.(2005)| . The Ice Cerenkov detector 
Ice Cube [Kestel et al.(2004)| at the South Pole and its 1 km 2 surface air 
shower detector Ice Top [Gaisser et al.(200 3)] are under construction. In Jan- 
uary 2008 40 Ice Cube strings and 40 surface detectors have been deployed, 
which implies the set-up is already 50% completed. Further experiments are 
the Telescope Array |Fukushima et al.(2007)| and its low energy extension 
TALE, as well as extensions of the Pierre Auger Observatory to lower en- 
ergies [Kalge s~et al.(2007)] . With this new high-quality data more detailed 
information will be available on the energy spectrum and the composition 
and it should be possible to distinguish between different scenarios for the 
transition from galactic to extragalactic cosmic rays. 

A promising complementary detection method for high-energy cosmic rays 
is the measurement of radio emission from air showers. This method provides 
three-dimensional information about air showers, similar to the fluorescence 
technique, but with the advantage of a much higher duty cycle. In the next 
years air showers are expected to be detected with the LOFAR radio obser- 
vatory |LOFA R(2007) |. An extensive research and development program is 
conducted in the Pierre Auger Collaboration with the goal to build a 20 km 2 
radio antennae array |van den Berg et al.(2007)| . 

The southern site of the Pierre Auger Observatory covers only a part of 
the whole sky, see Fig. [TBI Since the distribution of matter in the Universe 
is different in the parts seen from the northern and southern hemispheres 
it is important to observe the whole sky. The Northern Auger Observatory 
is designed to complete and extend the investigations begun in the South 
|Nitz et al. (2007)] . To unambiguously identify the sources of the highest- 
energy cosmic rays requires collecting many more events in spite of the steeply 
falling energy spectrum. The planned Northern site will be located in South- 
east Colorado, USA, having an instrumented area several times the area of 
Auger South. 

The Northern Observatory needs unrestricted support now, it is the next 
step in exploring the high-energy Universe in the upcoming years. With the 
completed Pierre Auger Observatory, with its Southern and Northern sites 
operated simultaneously, an exciting future in astroparticle physics is ahead 
of us. It will establish charged particle astronomy on the whole sky and will 
provide high-accuracy data to test astrophysical models of the origin of ultra 
high-energy cosmic rays. In addition, it will improve our understanding of fun- 
damental physics. The data will give insight into topics like the existence of 
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vacuum Cerenkov radiation, the smoothness of space, and tests of Lorentz in- 
variance | Klinkhammer and Risse(2008), Galaverni and Sigl(2007)|. Already 
now, the existing Auger data set stringent limits on theories. 
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